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Directed flow in ultrarelativistic nucleus-nucleus collisions is analyzed using the event plane from 
elliptic flow, which reduces the bias from nonflow effects. We combine this method with the de- 
termination of elliptic flow from Lee- Yang zeroes. The resulting method is more consistent and 
somewhat easier to implement than the previously used method based on three-particle cumulants, 
and is also less biased by nonflow correlations. Error terms from residual nonflow correlations are 
carefully estimated, as well as statistical errors. We discuss the application of the method at RHIC 
and LHC. 
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I. INTRODUCTION 



The standard method of analyzing directed flow (vi ) in nucleus- nucleus collisions using an estimate of the reaction 
plane [1, 2] was shown to be inadequate at ultrarelativistic energies [3], due to the smallness of directed flow and 
the relatively large magnitude of "nonflow" effects [2, 4]. Trivial nonflow effects such as momentum conservation can 
be taken into account in the analysis [3], and the standard analysis can be modified for this purpose [5]. A more 
systematic, model-independent way of eliminating nonflow effects was introduced [6], based on the observation that the 
three-particle average (cos(0i + 02 — 203)) (where 0i , 02 and 03 denote the azimuthal angles of three particles emitted 
in a collision, and (• • •) denotes an average over triplets of particles and events) is much less sensitive to nonflow effects 
• than the two-particle correlation (cos(0i — 02)) used in the standard analysis. Since (cos(0i + 02 — 203)) oc vfv2, 
where the elliptic flow V2 is large at ultrarelativistic energies, this provides an alternative way of analyzing directed 
flow. This method was first implemented at the CERN SPS by the NA49 Collaboration [7, 8], and recently led to 
the discovery of directed flow at RHIC [9]. Directed flow at ultrarelativistic energies is interesting in itself [10]; in 
addition, it is the only way of measuring the sign of elliptic flow, and to check experimentally that it is positive [11]. 

The practical implementation of the three-particle method presented in Ref. [6] is rather cumbersome and requires 
to estimate elliptic flow V2 independently, using cumulants [12] to avoid nonflow effects. In this paper, we suggest to 
analyze simultaneously elliptic and directed flows using the recently introduced method of Lee- Yang zeroes [13, 14] 
to analyze elliptic flow. This minimizes the bias from nonflow effects for both v\ and i>2- The practical recipe is 
presented in Sec. II. The theoretical background is briefly discussed in Sec. III. In Sec. IV, we derive the general 
order of magnitude of the systematic error due to nonflow correlations, which was underestimated in Ref. [6] . We also 
give analytical expressions of statistical errors. 



II. IMPLEMENTATION 



Let us first define useful quantities and notations. For a given event, we define the following complex-valued function 

M 

9l u92 {z) = II [T- + zewi(j)™s{<f> j -6 1 ) + zw 2 (j)cos(2(<i> j -62))] (1) 

where 9\ and 62 are angles, e is a real parameter, z is a complex variable, <pj are the azimuthal angles of the particles 
(measured using a fixed reference in the laboratory), and the product runs over all detected particles. As in other 
methods of analysis, w\ and W2 are weights appropriate to directed and elliptic flows, respectively, and can be any 
functions of particle type, transverse momentum pt and rapidity y. In Eq. (1), w n (j) is a shorthand for w n (pTj, y-j)- 
The best weight is the flow itself [12, 15], w n (pT, y) = v n (PT, y), where v n (pt, y) denotes the value of the flow in a small 
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(px, y) bin. In practice, one can choose as a first guess the center-of-mass rapidity for directed flow, w\ = y — ycM, 
and the transverse momentum for elliptic flow w 2 = Pti m regions of phase space covered by the detector acceptance. 
The average of g 01 ' 62 (z) over events in a centrality bin will be denoted by (g® 1 ' e2 (z)) cvts - The same notation holds for 
any observable associated with the event. 

The implementation of the method consists of three logical steps, as in the 3-particle correlation method [6]. The 
first one is to analyze integrated elliptic flow V 2 - By "integrated," we mean that it is summed over all detected 
particles with appropriate weights, and averaged over many events: 

V 2 = (y, Mi) cos(2(0 j - , (2) 

" = 1 'evts 

where the sum runs over all particles detected in an event and is the azimuthal angle of the reaction plane of the 
event. 1 The second step is the analysis of integrated directed flow V\, which is defined similarly: 

yi = /5>i(j)a>s(&-* fl )\ . (3) 

= 1 'evts 

Finally, in a third step, integrated values are used as a reference to analyze directed flow differentially, as a function 
of transverse momentum and rapidity. In practice, the method requires only two passes through the data: the first 
pass corresponds to Sec. II A, while the computations in Sec. II B and Sec. II C can be done in a single second pass. 



A. Integrated elliptic flow 

The analysis of integrated elliptic flow is identical to that presented in Refs. [13, 14], to which we refer the reader 
for further practical details. Note, however, that it uses the generating function introduced in Ref. [14], not that of 
Ref. [13]. Let us just give a short reminder of the recipe. 

One must evaluate g e t 1 ' 02 (z) for e = 6\ = 0, for several (typically 4 or 5) equally spaced values of 9 2 between and 
tt/2 (i.e., 9 2 = 0, 7r/(2n), . . . , (n — l)w/(2n) with n = 4 or 5) and many values of z on the imaginary axis, z = ir. One 
then plots the modulus \(g^' 62 {ir)) cvts \ for real, positive r as a function of r, for a fixed value of 9 2 . One determines 

0. Q 

numerically the value of r corresponding to the first minimum of this function, r = r 2 . The value of r 2 will be 
used in the analysis of directed flow below. It is closely related to elliptic flow, which can be estimated through the 
following formula [13, 14]: 

<4) 

where joi — 2.40483 is the first root of the Bessel function of the first kind Jo(x). The notation (• • -) e means an 
average over 9 2 . The analysis could in principle be performed with a single value of 6 2 . The averaging is only of 
practical importance: it reduces the statistical error on V 2 by a factor of ~ 2. 

Note that Eq. (4) yields an estimate of the absolute value of V 2 , not of V 2 itself. The sign of V 2 is determined 
simultaneously with directed flow, as we shall see below. 



B. Integrated directed flow 

The analysis of integrated directed flow involves the generating function g 01 ' 02 (z). The angle 9 2 takes the same 
lues as in Sec. II A. For a given 9 2 , z takes only one value, z = ir 02 . 

For each event, one first evaluates the generating function for e = 0\ = 0, as well as its derivative with respect to z: 



d 9o 02 ( -e 2] _ o.e 2( -e 2] ^ w 2 (j) cos(2(^- - 9 2 )) 
Oz T^i 1 + lr o w 2(j) cos(20j - 6 2 )) 



1 &r is the exact reaction plane, not an estimated plane; is unknown on an event-by-event basis. 



3 



One then evaluates g 
condition 



01,02 



Vg 2 ) for a single non-zero value of e, which is arbitrary but must satisfy the following 



V 2 

e<< W 



(6) 



Since V\ is a priori unknown, one must guess a reasonable value and check afterwards that the condition is satisfied. 
If e is too small, numerical errors may arise as the procedure amounts to expanding numerically g ei ' 02 (z) to order 
e 2 . Thus we recommend to perform tests with several values of e, and to check the stability of the results, before 
doing the full analysis. The angle d\ takes 5 (or more) equally spaced values of 9\ between and 27r, (i.e. Q\ = 
0, 2n/n 7 . . . , 2(n — l)n/n with n > 5). Note that the range differs from that of 82 values. 
Our estimate of V\ is defined by 



{V 1 fsgn(V 2 )= - 
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(cos(2(<?i -0 2 ))g^(ir e J)) 

\ I 0i, evts 



0,02 



V 



dz 



-C< 2 ) 



(7) 
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where e is an arbitrary small number satisfying condition (6), (• • -^.cvts denotes an average over 9\ and events, and 
Re is the real part of the ratio. Finally, sgn(V2) denotes the sign of V2, which is determined by this analysis of V\. 
The sign of V\, on the other hand, is not measured and must be postulated, as in any other method of analysis. 

As in Sec. II A, the averaging over 62 in Eqs. (7) and (9) is only for practical purposes: it reduces statistical errors 
by a factor of <~ 2. In contrast, the averaging over 6\ cannot be avoided, for reasons we shall explain below in Sec. III. 

C. Differential directed flow 

One can then turn to the analysis of differential flow, i.e., the flow of particles of a given type in a definite phase- 
space window, which we shall call "protons" for the sake of brevity. A "proton" azimuth will be denoted by ip, and 
the corresponding differential directed flow by v[. 

The estimate of v[ involves the derivative of the generating function with respect to the proton weight, evaluated 
at z = irn 2 : 



dd! 1 ' 92 ^ 2 ) = 1 ,02 (ir 02 ) ir^cos^-flQ 

dwity) ye {o l + ir e Q 2 ewi(V>)cos(V> - 0i) + ir e 2 w 2 {ip) cos(2(i{> - 9 2 ))' 



(8) 



and is defined by 



dw\(ip) 



i,4> 



0,6» 2 



V 



dz 



M 2 ) 



(9) 



where (• • ■)e 1 ^ denotes an average over protons and over values of Q\. Note that the value of integrated directed flow 
V\ is only required after averaging over protons or events, which means that the ratio between parentheses can be 
computed at the same time as that in Eq. (7). Moreover, the denominator of this ratio is the same as in Eq. (7), so that 
in the second pass through data, one only need compute three quantities for each event: this common denominator, 
the numerator of Eq. (7) (actually, a value for each angle 6*i), and the numerator of Eq. (9), taking into account only 
the "protons." 



D. Relation with the three-particle method of Ref. [6] 



The three-particle method we proposed in Ref. [6] is based on averages of the type (cos(</>i + </> 2 — 2</> 3 )), where 
4>i, <f>2 and 03 are the azimuthal angles of three particles belonging to the same event. Let us explain how such 
three-particle averages appear in the present method: expanding g^ 1 '° 2 (z) in Eq. (1) to order e 2 z 3 produces terms of 
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the type cos(</>i — 9\) cos(</>2 — 9i) cos(2(</>3 — 02))- Multiplying by cos(2(6*i — 6*2)) and averaging over 9\ and 62, as in 
Eqs. (7) and (9), one obtains | 008(0! + 2 — 2</> 3 ), thus recovering the three-particle averages. The present method 
uses an expansion to order e 2 (see Sec. Ill), but for a fixed value of z, namely ir® 2 : it therefore also includes the 
information from higher-order correlations, generated by higher powers of z, not only three-particle correlations. 



III. THEORY 



The general philosophy is that zeroes of generating functions are direct probes of collective effects [13], while they 
are little sensitive to nonflow effects. The key property of the generating function in Eq. (1) is its factorization 
property: if a nucleus- nucleus collision can be viewed as the superposition of independent subsystems containing a 
few particles (for instance, independent nucleon-nucleon collisions), i.e., if there are no collective effects, the average 
over events (g® 1 ' 02 (z)) cvts is the product of the contributions of individual subsystems. The zeroes of the generating 
function are then simply the zeroes of the contributions of the subsystems. Since the value of the generating function 
is 1 for z = by construction, the value of z for which it vanishes verifies \z\ ~ 1 for dimensional reasons (in this 
discussion, we assume unit weights w\ = W2 = 1 for simplicity). 

If, on the other hand, there is collective flow, nothing prevents the zeroes from moving much closer to the origin, 
\z\ <C 1, and this is indeed what happens: the positions of the zeroes scale with the event multiplicity like 1/M (while 
they would be independent of M for independent subsystems). This relates our method to the Lee- Yang theory of 
phase transitions [16], where the zeroes of the grand partition function move closer and closer to the real axis as the 
size of the system increases when, and only when, there is a phase transition. 

Let us now be more quantitative. We first evaluate the average value of g e t 1 ' e2 {z) for events having exactly the same 
reaction-plane orientation <fr R (denoting such an average by (• • • |$r)). Taking the logarithm of Eq. (1) and expanding 
to order z, we obtain 

In (g e ^ 2 {z)\<S> R ) ~ ze ^ Wl (j) cos(& - + z (j^ Mj) cos(2(^ - 2 ))|<^ 

~ zeVicos($ R -6 1 ) + zV 2 cos(2(<f> R -62)), (10) 

where we have used Eqs. (2) and (3) and assumed symmetry with respect to the reaction plane, which implies 
(12 j w z(j) sin(2(0j — $fi)))evts = (12 j w i(j) sin((j)j — <&ft)) ev ts = 0. Terms of order z 2 and higher in the expansion arc 
responsible for systematic errors that will be estimated in Sec. IV A. The average over events is eventually obtained 
by averaging over which is randomly distributed: 

W*W)^ = f*W*i*)\**)^- (ID 

For e = 0, substituting the estimate (10) into Eq. (11), one obtains 

[9o' e2 (z)) =Io(\V 2 \z), (12) 

\ / evts 

where I is the modified Bessel function of order 0. Its zeroes lie on the imaginary axis, the first one being located at 
ijoi- This explains how | V2 1 is obtained in Sec. II A, Eq. (4). Since zeroes are expected to be on the imaginary axis 
except for irrelevant statistical fluctuations and/or detector effects [13], we suggest to look for the position of the first 
minimum of K^g'^ 2 (ir)) cvts | on the upper imaginary axis, denoted by z = ir® 2 in Sec. II A, rather than that of the first 
zero of (g^'® 2 (z)) cvts in the complex plane. 

In order to obtain an estimate of V\, we shall now study how the first zero of (gf 1 ' e2 (z)) cvts , which we denote by 
z e , varies for small e. For small enough e, we may write 

_ {g e y 2 (zo)) cvts 

where zo is the value of z t for e = 0. According to our general philosophy, the value of z e is directly related to 
collective flow, and little sensitive to nonflow effects. Then, this is also true of the right-hand side (r.h.s.) of Eq. (13), 
which is experimentally measurable. Following the discussion above, we replace the exact zero z by irQ 2 . To relate 
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the quantity in Eq. (13) to anisotropic flow, we use our theoretical estimate of the generating function, defined by 
Eqs. (10) and (11). The denominator is simply obtained by differentiating Eq. (12): 

\^~ (Zo) ) =N/i(IW. (14) 

\ / evts 

The numerator is obtained by substituting Eq. (10) into (11), and expanding to order e 2 : 

<S^o)) cvts = j4h(V 2 z ) V? 003(2(6, -6 2 )). (15) 

Since I\ is an odd function, we may write I\(V 2 z ) = V^l^o) sgn(V 2 ). Putting together the last two equations, 
replacing z with ir-g 2 and | V2 1 with its estimate joi/ r o 2 (see Eq. (4)), we obtain 

V?sgn{V 2 )cos{2{e 1 -e 2 )) = -^[r e A (16) 

This yields in principle an estimate of V\ for each value of Q\ and 6 2 . However, the expansion of (g e e 1,02 (ir^ i )) ev ts in 
powers of e generally yields a non-vanishing term of order e due to statistical fluctuations. Multiplying the previous 
equation by cos(2(#i — 2 )), and averaging over 5 or more equally spaced values of 6\, one eliminates this term. 
Averaging over 9 2 , one obtains our final estimate, Eq. (7). 

We can now turn to differential directed flow. For this purpose, let us study how the zero of g^ 1 ' 02 (z) moves when 
the weights w\ of all protons are shifted by some small quantity 5w\ : 

dg e y 2 (z*) 



dwi{ip) 



z Sw -z = — 8, ' ~r~ Swi - ( 17 ) 

-(zo) 



z 



The value of z$ w is directly related to collective effects. Following the same procedure as for integrated directed flow, 
we estimate the value of the r.h.s. when flow is present. To evaluate the numerator, we first compute the average for 
a fixed Using Eq. (8), assuming that the proton and the other particles are uncorrelated for fixed (i.e., that 
the correlation between the proton and the other particles is only due to flow), and neglecting the contribution of the 
proton to the generating function, we obtain 



\$r^ - eMcos(V>-0i)|$fl)(<7^ 2 (zo)|$ii) 



dwi(tp) 

= ez v' 1 cos($R-Oi)exp[zoeVicos($R-6i)+zoV2Cos(2($ R -6 2 ))], (18) 

where we have used Eq. (10), the definition of v[ = (cos^ — $r)), and assumed symmetry with respect to the reaction 
plane, (sin(^ — $_r)) = 0. We then expand the exponential to order e and integrate over $_r: 

(), *<"" 2{Z{,] \ - e2 z * v > lVl /i(V 2 « )cos(2((9i-fl 2 )). (19) 



d Wl {ip) I 2 
Using Eq. (14), and replacing z with ir-g 2 , we obtain 



v' lVl sgn(U 2 ) cos(2(e 1 - 6 2 )) = ) ^ W A . (20) 

\ / cvts 

Multiplying both sides of this equation by cos(2(#i — 9 2 )) and averaging over 61, one eventually obtains Eq. (9). 

For the sake of consistency, one must recover the integrated flow V\ by integrating the differential flow v[ over phase 
space. Our estimates (7) and (9) satisfy this sum rule, provided e is small enough. In order to prove this, we weight 
Eq. (9) with w\(ip) and integrate over phase space. The following quantity appears 



M 
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Next, we use the fact that g 01 ' 02 ^ 02 ) scales like e 2 for small e, hence 




(22) 



This completes the proof. 



IV. ERRORS 



The main motivation for analyzing collective flow with multiparticle correlations is that one thereby reduces spuri- 
ous, "nonflow" effects. The magnitude of errors due to residual nonflow effects is estimated in Sec. IV A. The price to 
pay for this greater reliability is an increase in statistical errors, which are evaluated in Sec. IV B. Finally, azimuthal 
asymmetries in the detector acceptance are a potential source of bias in the analysis. They are carefully studied in 



There are several tricks to check experimentally whether or not the analysis is biased by nonflow effects. A first 
one is to perform the analysis twice, using different weights in Eq. (1), e.g., one analysis with W2 — pt an d another 
with w 2 = 1, or redoing the analysis with zero weights for the particles in one hemisphere. The final differential flow 
results should be independent of the weights, while no such independence is expected for nonflow effects. However, 
using non-optimal weights increases statistical errors, sometimes by large amounts, so that this is doable in practice 
only if statistical errors are not a limitation. A second trick is to check that the final results look "reasonable:" this 
means in particular that they must be consistent with the known symmetries of the system, i.e., v\{—y) = —vi(y), 
v 2(—y) = V2(y) for a symmetric collision. Nonflow effects can destroy these symmetries: for instance, effects of 
total momentum conservation may contaminate the measurement of v±, yielding an estimate that does not vanish 
near midrapidity [5]. This may however occur only when the weights in Eq. (1) are not symmetric themselves (i.e., 
w i(—y) 7^ ~ w i(y), w 2{—y) ^ w 2(y))- With symmetric weights, the fact that v\ vanishes at midrapidity is not an 
indication that it is not biased by nonflow effects. 

In this paper, we are interested in directed flow at ultrarelativistic energies, which is most often at the border of 
observability, so that it is not always possible to use the above tricks. It is therefore important to derive estimates of 
the magnitude of nonflow effects using purely theoretical arguments. There is of course no way to derive quantitative 
estimates, due to the variety of the physical effects involved [3]. In most cases, one can at best rely on orders 
of magnitude and simple scaling rules [12], which are derived below. More quantitative statements can be made 
concerning momentum conservation, a well-known bias in analyses of directed flow [5, 17, 18]. 

With standard methods of flow analysis [1, 2] and also, to a lesser degree, with cumulants [12], the analysis may 
yield a non-zero value of v\ and t>2 even when they are in fact zero, i.e., even if only "nonflow" effects are present [19]. 
This is not the case with Lee- Yang zeroes, where a non-vanishing result beyond statistical errors can be considered 
a clear signal of anisotropic flow. 

However, the interference of collective flow and nonflow effects may produce a small relative error, which is the sum 
of two terms: 



Let us explain how such terms arise. Our determination of directed flow with the present method involves three- 
particle averages such as (cos(0i + <f> 2 — 2</> 3 )) (see Sec. IID). 2 The contribution of flow to this average is v\v 2 - Let us 
now evaluate the contribution of nonflow effects. As a simple model, assume that the M particles detected in an event 
are emitted in M/2 collinear pairs ("jet-like correlations"). There is a probability 1/(M — 1) ~ l/M that cf>2 — <f>i, in 
which case the three-particle average becomes (cos2((/>i — fa)) = v\. This gives the first term in Eq. (23). Similarly, 
there is a probability of order l/M that <\>\ — 03, in which case the three-particle average reads (cos(</>i — fa)) = vf; 
this yields the second term in Eq. (23). Since the present method is expected to be useful essentially when v\ < 



Sec. IV C. 



A. Systematic errors from nonflow effects 




(23) 



2 One can check that the systematic errors due to higher-order correlations arc at most of the same order as those due to three-particle 
correlations. 
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this term is subleading. Finally, one can check explicitly that momentum conservation does not contribute to the 
error terms in Eq. (23). 

The order of magnitude of the systematic error, Eq. (23), is more general than suggested by this toy model of 
nonflow correlations. It can be proven more rigorously by pushing the expansion of the logarithm of the generating 
function, Eq. (10), to order z 2 , and studying the influence of the additional terms on the flow estimate, much in 
the same way as in Ref. [13], where the interference between flow and nonflow effects was also considered (see the 
discussion following Eq. (48) in the reference). Of course, Eq. (23) is but a scaling law, and each term involves an 
unknown numerical coefficient. One can try to estimate these coefficients experimentally by studying the difference 
between flow estimates from various methods, using situations where this difference is believed to be dominated by 
nonflow effects. This strategy was applied by the STAR Collaboration [9] in their estimates of systematic errors on 
w 4 and v\. 

With the standard, event-plane method [1, 2], the error due to nonflow effects is much larger: 

(s VlUnRow=o( i y 



V! \Mv{ 

Comparing with Eq. (23), our method reduces the nonflow error by a factor of l/t>2, i.e., typically 15 to 20 at RHIC. 

The errors with the present method, Eq. (23), also apply to the three-particle cumulant method, but they were 
not mentioned in Ref. [6]. The cumulant method has an additional error term, coming from the "pure nonflow" 
contribution to the three-particle average (cos(</>i + (p 2 — 2^3)), of order 1/M 2 . This gives an additional error term in 
Eq. (23) 

(5«l)nonflow _ Q ( 1 Y (2g) 



Vl \M 2 V2Vi 

This contribution, which was the only error term mentioned in Ref. [6] , disappears with Lee- Yang zeroes. In practice 
at ultrarelativistic energies Mv 2 > 1, so that the additional term is typically of the same order of magnitude as the 
first term of Eq. (23). In addition, momentum conservation produces an error of this order, so that the Lee- Yang 
zeroes method does a better job in eliminating effects of momentum conservation. Finally, the three-particle method 
requires to estimate v 2 independently. If v 2 is estimated from two-particle cumulants or by the event-plane method, 
the resulting error on v\ is (<5vi) non flow/wi = 0(1/Mv 2 ), which may dominate over both terms in Eq. (23). 



B. Statistical errors 



Statistical uncertainties on V\ and v[ naturally involve both directed and elliptic flows. Both errors turn out to scale 
like the statistical error on the differential elliptic flow v' 2 determined from Lee- Yang zeroes [13], which is natural since 
we use elliptic flow as a reference. The statistical error on elliptic flow involves the associated resolution parameter 

V? 

X2 = -, (26) 

where er 2 is a measure of event-by-event fluctuations. Provided nonflow effects are not too large, (72 can be obtained 
experimentally by the following formula: 

I m \ 

(27) 

More accurate determinations are discussed in Ref. [13]. Once \i nas been determined, the statistical error on v' 2 
is [13]: 

(5v' 2 f = 4JVV1U1) 2 ( 6XP (11 ° OS202 ) J o(2.7oisin0 2 ) - exp ^-||cos20 2 ) J (2j 01 cos0 2 )^ , (28) 

where N' is the number of "protons" in the phase-space region under study. For large \2, this yields Sv' 2 = 1/V2-/V', 
which is the expected statistical error when the reaction plane is exactly known. When \i becomes significantly 
smaller than unity, on the other hand, the statistical error on v' 2 increases exponentially as Sv' 2 oc exp(2.9/x 2 ), an d 
the method cannot be applied [13]. 
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Let us now discuss statistical errors on the directed flow estimates. These errors involve two parameters o\ and \i 
jfined as in Eqs. (26) and (27), ] 
ssociated with directed flow. 3 In 
following substitution in Eq. (28): 



defined as in Eqs. (26) and (27), namely, a\ ~ (VV w|) ev ts and Xi — Vi/&i- The latter is the resolution parameter 
associated with directed flow. 3 In order to obtain the square of the statistical error on V^sgn^), one then does the 



J_ ^ 2,1(1 + X j) 

N> N cvts ( ' 

For large xi and X2> this yields 8V\ = a\/^2N evts , which is the expected value of the statistical error if the reaction 
plane is exactly known. 

Finally, the statistical error on the differential directed flow v[ is simply related to the error on v' 2 : 

(H) 2 = 1 -^{5v' 2 )\ (30) 

Al 

One notes that the statistical error on directed flow obtained with this method is always larger than the statistical 
error on elliptic flow in absolute value. When \2 is large enough, Sv^ is the same as with the event-plane analysis. 

When using 3-particle cumulants [6, 9], the statistical error is given by a similar formula, where 5v' 2 is the error on 
elliptic flow obtained with the event-plane method (or by 2-particle correlations). Going from the 3-particle method 
to the present method, the increase in statistical errors will be the same as when going from the standard method to 
Lee- Yang zeroes in the analysis of elliptic flow. For a semi-central Au-Au collision at RHIC analyzed with the STAR 
detector, the resolution parameter \2 is at least 1, which means a factor of at most 2 increase in statistical errors [13]. 
This increase is compensated by the smaller error from nonflow effects. 



C. Azimuthal asymmetries in the detector acceptance 

A nice feature of all flow analyses based on cumulants or Lee- Yang zeroes (which amount to cumulants of many- 
particle correlations [13]) is that they automatically eliminate most effects of azimuthal asymmetries in the detector 
acceptance: more precisely, such asymmetries cannot produce a signal by themselves, even if the detector has a very 
partial azimuthal coverage. This is not the case with the event-plane analysis, where several flattening procedures 
are required [2], and with two-particle correlation methods [20, 21], which use mixed events to correct for acceptance 
effects. 

Acceptance asymmetries generally have two effects with cumulants or Lee- Yang zeroes [12, 13]: 1) the flow given 
by the analysis differs from the true flow by some factor, which can be computed analytically once the acceptance 
profile is known; 2) different harmonics interfere, so that a measurement of v\ is generally biased by v 2 and vice- versa. 
These interference terms can also be calculated analytically. 

Such interference terms between v\ and v 2 are also present here, and they turn out to be a rather serious problem. 
Indeed, as an explicit calculation will show, acceptance asymmetries produce an additional term proportional to v 2 
in the left-hand side (l.h.s.) of Eq. (7): even if there is only elliptic flow in the system, the analysis then yields a 
spurious directed flow v\ cx ^/v^. Since we are typically interested in a situation where v 2 is large and v\ is small, 
such an interference is most unwelcome, and it must be thoroughly studied. 

In order to focus on the interference with v 2 , we assume that only elliptic flow is present in the system, Vi = 0, and 
we evaluate the spurious directed flow. The following study is limited to the integrated flow V\. For a given reaction 
plane $_r, the azimuthal distribution of a particle of type j is, before detection, 

dN 

— cx 1 + 2v 2 (j) cos(2(& - $«)). (31) 
With this distribution, and with a non-symmetric detector, the first term in Eq. (10) becomes 

(j2 cos( ^' ~ = Re [^o acc e^ ei + V r 2 acc e i(2 * H - ei) + y 2 '-c c i(-2* J? -e l) ] ) (32 ) 



3 Please note that the resolution parameter \ as defined in Ref. [2] is larger by a factor of \/2 
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where Vq cc , U 2 acc and V 2 ' acc are complex coefficients defined by 



M 



v o cc = (E^o> i0 y 

/ M 



\j=i i 

V^ cc = (jrw^v^c^y (33) 

with the averages taken over detected particles. For a symmetric acceptance, the three coemcients vanish. With 
these values, one can compute the generating function (gf 1 ' e2 (z)) to order e 2 following essentially the same steps as 
in Sec. 111. Let us give the final result, skipping tedious intermediate calculations: to leading order in acceptance 
asymmetries, one must replace (Ul) 2 with Re pV^f^Kf 00 ] in the l.h.s. of Eq. (7). 

For moderate inhomogeneities in the detector acceptance, the interference of elliptic flow on the estimate of directed 
flow remains under control. Thus, assuming an elliptic flow value v-i — 0.06, variations of 20% in the detector efficiency 
will result in a spurious v\ of at most 0.01. Using appropriate weights (e.g., weighting with the inverse of the detector 
efficiency profile) would significantly decrease this spurious directed flow value. However, we fear that working with 
a detector with incomplete azimuthal coverage might prove prohibitive. 

Finally, note that the three-particle method is affected by acceptance asymmetries in essentially the same way, but 
the corresponding interference term was unfortunately omitted in Ref. [6] . 



V. CONCLUSIONS AND PERSPECTIVES 



In the foregoing, we have introduced a new method to analyze directed flow v\ in cases where it is too small to be 
reliably obtained with standard two-particle methods without contamination from nonflow effects, and also too small 
to be determined by more recent methods as cumulants or Lee- Yang zeroes, which would both give large statistical 
uncertainties. Moreover, the method relies on the implicit assumption that elliptic flow is reasonably large, as when 
using three-particle cumulants [6]. With respect to the latter, the new method is easier to implement: instead of 
interpolating successive derivatives of a function, one need just determine the position tq of the first minimum of a 
function, and then compute the values at ro of the three quantities that appear in the numerators and denominator 
of Eqs. (7) and (9). In addition, the method is conceptually more elegant, since it relies on the deep relation there 
exists between the behaviour of the zeroes of generating functions and the presence of collective effects in the systems 
these functions describe [13]. 

We have also carefully estimated both systematic errors due to nonflow effects and statistical uncertainties arising 
from finite available statistics. We may use these estimates to discuss the applicability of the method to present and 
future heavy-ion experiments. At RHIC, the three-particle cumulant method was already successfully applied [9]. In 
mid-central collisions, statistical uncertainties should be at most a factor 2 larger with the new method, but with the 
high statistics of Run 4, this should not be a problem. On the other hand, systematic errors will be reduced: the error 
term (25), which arises in particular from momentum conservation, disappears with the new method. The method 
could be useful for new measurements of directed flow at the CERN SPS, although this would require detectors with 
a larger coverage than what single experiments had in the past runs. 

Finally, one can expect that the present method will allow measurements of directed flow at LHC, if any, using 
muons from decay pions seen in the ALICE spectrometer at forward rapidities (2.5 < rj < 4) or hits in the CMS 
very forward hadronic calorimeter (3 < rj < 5). Since V2 at LHC is expected to be at least as large as at RHIC, 
and the multiplicity will be higher, elliptic flow will be analyzed with an excellent resolution. Then, the statistical 
errors with our new method will be barely larger than with the standard event-plane method: statistics will not be 
a limitation down to values of v\ of a fraction of a percent. Systematic errors from nonflow effects may be a more 
severe problem, although our method minimizes their magnitude. Comparing Eq. (30), with xi ~ Vi\/~M, and the 
first term of Eq. (23), one sees that as v\ decreases, systematic errors tend to become larger than statistical errors. 
One may reasonably hope, however, that the huge particle multiplicity M expected at LHC will compensate for the 
smaller value of V\, and that V\ will eventually be observed. 
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